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A general method exists for studying Abelian and non-Abelian gauge theories, as 
well as Euclidean quantum gravity, at one-loop level on manifolds with boundary. 
In the latter case, boundary conditions on metric perturbations h can be chosen to 
be completely invariant under infinitesimal diffcomorphisms, to preserve the invari- 
ance group of the theory and BRST symmetry. In the de Dondcr gauge, however, 
the resulting boundary-value problem for the Laplace type operator acting on h is 
known to be self-adjoint but not strongly elliptic. The present paper shows that, 
on the Euclidean four-ball, only the scalar part of perturbative modes for quan- 
tum gravity is affected by the lack of strong cllipticity. Interestingly, three sectors 
of the scalar-perturbation problem remain elliptic, while lack of strong cllipticity 
is "confined" to the remaining fourth sector. The integral representation of the 
resulting ^-function asymptotics on the Euclidean four-ball is also obtained; this 
remains regular at the origin by virtue of a peculiar spectral identity obtained by 
the authors. There is therefore encouraging evidence in favour of the £(0) value 
with fully diff-invariant boundary conditions remaining well defined, at least on 
the four-ball, although severe technical obstructions remain in general. 
2000 MSC. Primary 58J35, 83C45; Secondary 81S40, 81T20. 



1. Introduction 

This paper is motivated by the authors' struggle over many years with 
an important problem in quantum field theory and spectral geometry, i.e. 
the functional determinant in Euclidean quantum gravity on manifolds with 
non-empty boundary. The related open issues are not yet settled, but there 
is a sufficient amount of new calculations to justify further efforts, as we 
are going to see shortly. 

The subject of boundary effects in quantum field theory (Deutsch and 
Candelas [1]) has always received a careful consideration in the literature 
by virtue of very important physical and mathematical motivations, that 
can be summarized as follows. 

(i) Boundary data play a crucial role in the functional-integral approach 
(DeWitt [2]), in the quantum theory of the early universe (Hartle and 
Hawking, Hawking [3]) in supergravity (Hawking [4]) and even in string 
theory (Abouelsaood et al. [5]). 

(ii) The way in which quantum fields react to the presence of boundaries 
is responsible for remarkable physical effects, e.g. the attractive Casimir 
force among perfectly conducting parallel plates (Bordag et al., Milton, 
Nesterenko et al. [6]), which can be viewed as arising from differences of 
zero-point energies of the quantized electromagnetic field. 

(iii) The spectral geometry of a Ricmannian manifold (Gilkcy [7]) with 
boundary is a fascinating problem where many new results have been de- 
rived over the last few years (Kirsten [8], Vassilevich [9]). 
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(iv) Boundary terms (Moss [10]) in heat-kernel expansions have become 
a major subject of investigation in quantum gravity (Avramidi [11]), since 
they shed new light on one-loop conformal anomalies (Esposito et al., Moss 
and Poletti [12], Tsoupros [13]) and one-loop divergences (Esposito [14], 
Esposito ct al. [15]). 

In our paper we are interested in boundary conditions for metric pertur- 
bations that are completely invariant under infinitesimal diffeomorphisms, 
since they are part of the general scheme according to which the boundary 
conditions are preserved under the action of the symmetry group of the 
theory (Barvinsky [16], Moss and Silva [17], Avramidi and Esposito [18]). 
In field-theoretical language, this means setting to zero at the boundary 
that part ttA of the gauge field A that lives on the boundary B (n being a 
projection operator): 



(1) 



ttA 

as well as the gauge-fixing functional, 

'*{A)] b = 0, (2) 

and the whole ghost field 

Mb = o. (3) 

For Euclidean quantum gravity, Eq. (1) reads as 

[hijh = 0, (4) 

where hij are perturbations of the induced three-metric. To arrive at the 
gravitational counterpart of Eqs. (2) and (3), note first that, under in- 
finitesimal diffeomorphisms, metric perturbations h^ v transform according 
to 

h^ = h^ + V^(p v ), (5) 

where V is the Levi-Civita connection on the background four-geometry 
with metric g, and <f u dx u is the ghost one-form (strictly, our presentation is 
simplified: there are two independent ghost fields obeying Fermi statistics, 
and we will eventually multiply by —2 the effect of tp v to take this into 
account). In geometric language, the infinitesimal variation <5/i M „ = h^ v — 
h^u is given by the Lie derivative along ip of the four- metric g. For manifolds 
with boundary, Eq. (5) implies that (Esposito et al. [19], Avramidi et al. 
[20]) 

= + ip^j) + Kijip , (6) 
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where the stroke denotes three-dimensional covariant differentiation tangen- 
tially with respect to the intrinsic Levi-Civita connection of the boundary, 
while Kij is the extrinsic-curvature tensor of the boundary. Of course, ifo 
and ipi are the normal and tangential components of the ghost, respectively. 
By virtue of Eq. (6), the boundary conditions (4) are "gauge invariant", 
i.e. 



hij 



(7) 



if and only if the whole ghost field obeys homogeneous Dirichlet conditions, 
so that 

bole = 0, (8) 

[<Pi]B = 0. (9) 

The conditions (8) and (9) are necessary and sufficient since ipo and ipi 
are independent, and three-dimensional covariant differentiation commutes 
with the operation of restriction to the boundary. We are indeed assuming 
that the boundary B is smooth and not totally geodesic, i.e. ^ 0. 

However, for totally geodesic boundaries, having Kij = 0, the condition (8) 
is no longer necessary. 

On imposing boundary conditions on the remaining set of metric per- 
turbations, the key point is to make sure that the invariance of such bound- 
ary conditions under the infinitesimal transformations (5) is again guaran- 
teed by (8) and (9), since otherwise one would obtain incompatible sets 
of boundary conditions on the ghost field. Indeed, on using the DeWitt- 
Faddeev-Popov formalism for the (out | in) amplitudes of quantum gravity, 
it is necessary to use a gauge-averaging term in the Euclidean action, of 
the form 2 

where $ u is any functional which leads to self-adjoint (elliptic) operators 
on metric and ghost perturbations. One then finds that 

s%(h) = %(h)-%(h) = r; <p v , (ii) 

where is an elliptic operator that acts linearly on the ghost field. Thus, 
if one imposes the boundary conditions 

* M W] B = 0, (12) 
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and if one assumes that the ghost field can be expanded in a complete 
orthonormal set of eigenfunctions of T£ which vanish at the boundary, 
i.e. 



(13) 



(A) 



(14) 



(A) 



o, 



(15) 



the boundary conditions (12) are automatically gauge-invariant under the 
Dirichlct conditions (8) and (9) on the ghost. 

Having obtained the general recipe expressed by Eqs. (4) and (12), we 
can recall what they imply on the Euclidean four-ball. This background 
is relevant for one-loop quantum cosmology in the limit of small three- 
geometry on the one hand (Schlcich [21]), and for spectral geometry and 
spectral asymptotics on the other hand [8, 9]. As shown in Ref. 19, if one 
chooses the de Donder gauge-fixing functional 

1 



(16) 



which has the virtue of leading to an operator of Laplace type on in the 
one-loop functional integral, Eq. (12) yields the mixed boundary conditions 



dh 



oo 



dr 



-h li 

.2 0j 



= 0, 



(17) 



dh M + 3^ _ 1 dhpo 
dr t 1 2 dx l 



= 0. 



(18) 



In Refs. 15, 19, the boundary conditions (4), (17) and (18) were used to 
evaluate the full one-loop divergence of quantized general relativity on the 
Euclidean four-ball, including all h^ v and all ghost modes. However, the 
meaning of such a calculation became unclear after the discovery in Ref. 
18 that the boundary-value problem for the Laplacian P acting on metric 
perturbations is not strongly elliptic by virtue of tangential derivatives in 
the boundary conditions (17) and (18). Moreover, the work by Dowker and 
Kirsten [22] had proved even earlier, in a simpler case, that the boundary- 
value problem with tangential derivatives is, in general, not strongly elliptic. 
Strong cllipticity [8, 18] is a technical requirement ensuring that a unique 
smooth solution of the boundary-value problem exists which vanishes at 
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infinite geodesic distance from the boundary. If it is fulfilled, this ensures 
that the L 2 trace of the heat semigroup e~ tp exists, with the associated 
global heat-kernel asymptotics that yields one-loop divergence and one-loop 
effective action. However, when strong ellipticity does not hold, the L 2 trace 
of e~ tp acquires a singular part [18] and hence (-function calculations may 
become ill-defined. 

All of this has motivated our analysis, which therefore derives in Sec. 
2 the eigenvalue conditions for scalar modes. Section 3 obtains the first 
pair of resulting scalar-mode (-functions and Sec. 4 studies the remaining 
elliptic and non-elliptic parts of spectral asymptotics. Results and open 
problems are described in Sec. 5. 



2. Eigenvalue conditions for scalar modes on the four-ball 

On the Euclidean four-ball, which can be viewed as the portion of flat 
Euclidean four-space bounded by a three-sphere of radius q, metric per- 
turbations hfjtv can be expanded in terms of hypersphcrical harmonics as 
(Lifshitz and Khalatnikov [23], Esposito et al. [24]) 



h 00 (x,T) = J2a n (T)Q^(x), 



(19) 



71=1 



h i(x,r) = 



n=2 



(X) 



b n {T)^^- )+ c n {r)S\ n \x 



(20) 



71—3 



g[g(*) 
(n 2 - 1) 



Ee„(r) 

n— 1 



{n) 



(X) 



n— 3 



[U(t)(s$(x) + S$(x)) + k n {r)G^{x)] , (21) 



where r G [0, q] and Q in) {x), S^' {x) and G\j>(x) are scalar, transverse vec- 
tor and transverse-traceless tensor hypersphcrical harmonics, respectively, 
on a unit three-sphere with metric c^. By insertion of the expansions 
(19)-(21) into the eigenvalue equation for the Laplacian acting on h^, and 
by setting \f~E — > iM, which corresponds to a rotation of contour in the 
(■-function analysis (Barvinsky et al. [25]) one finds the modes as linear 
combinations of modified Bcsscl functions of first kind. Modified Bessel 
functions of the second kind are not included to ensure regularity at the 
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origin r = 0. For details, we refer the reader to the work by Esposito et al. 
[26]. 

The boundary conditions (4), (17), (18), (8), (9), jointly with the mode- 
expansions on the four-ball, can be used to obtain homogeneous linear sys- 
tems that yield, implicitly, the eigenvalues of our problem. The conditions 
for finding non-trivial solutions of such linear systems are given by the van- 
ishing of the associated determinants; these yield the eigenvalue conditions 
8(E) = 0, i.e. the equations obeyed by the eigenvalues by virtue of the 
boundary conditions. For the purpose of a rigorous analysis, we need the 
full expression of such eigenvalue conditions for each set of coupled modes. 
Upon setting y/~E — > iM, we denote by D(Mq) the counterpart of 8(E), 
bearing in mind that, strictly, only 8(E) yields implicitly the eigenvalues, 
while D(Mq) is more convenient for ^-function calculations [25]. 

In particular, we here focus on scalar modes (for the whole set of modes, 
see again the work in Ref. 26). For all n > 3, coupled scalar modes 
a„,()„,(i„, e n are ruled by a determinant reading as 

D n (Mq) = det Pij (Mq), (22) 

with degeneracy n 2 , where pij is a 4 x 4 matrix with entries (hereafter, /„ 
are modified Bessel functions of first kind) 

p n = In(Mq) - MqI' n (Mq), p 12 = MqI' n (Mq), (23) 

p 13 = (2 - n)I n - 2 (Mq) + MqI' n _ 2 (Mq), 

p 14 = (2 + n)I n+2 (Mq) + M qI' n+2 (Mq) , (24) 

P2i = -(n 2 - l)In(Mq), p 22 = 2MqI' n (Mq) + 6I n (Mq), (25) 

p 23 = 2(n + l)MqI' n _ 2 (Mq) - (n 2 - 6n - 7)I n . 2 (Mq), (26) 

p 24 = -2(n - l)MqI' n+2 (Mq) - (n 2 +6n- 7)I n+2 (Mq), (27) 

P31 = 0, P32 = -In(Mq), (28) 

P33 = J^ ± Y) In - 2{Mq) ' P3i = J^TY) In+2{Mq) ' (29) 

p 4 i = 3I n (Mq), p 42 = -2I n (Mq), p 43 = -I n _ 2 (Mq), p 44 = -I n+2 (Mq). 

(30) 

The hardest part of our analysis is the investigation of the equation obtained 
by setting to zero the determinant (22). For this purpose, we first exploit 
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the recurrence relations among I n ,I n+ i and I' n to find (from now on, w 
Mq) 



pn = I n {w) - wl' n (w), pi2 = wl' n (w), pi 3 = wl' n (w) + nl n (w), 
Pi4 = wl' n {w) - nl n (w), 

P2i = -(n 2 - l)In(w), p 22 = 2(wl' n (w) + 3I n (w)), 

2n(n- l)(n + 3)" 



P23 = (n + 1) 



+ 2 



w + 



3(n + l) + 
(n- l)(n + 3) 



W 



P24 



= (n-l){ 



w + 



3(n-l) + 
(n + l)(n-3) 



2n(n + l)(n-3) 



/«(w) 



P33 = 



PU 



(n + 1) 
(n-2) 

("-!) 
(n + 2) 



P43 



P31 = 0, p 3 2 = -In(w), 

2n(n - 1) 



1 



2n(n + l) 



P4i = 3J„(u;), p 4 2 = -2I n (w), 

1 + 



2n(n-l)\ 2(n-l) 

—-5 — J i„H — i n (w), 



( 2n{n+l) \ 2(n+l) 



(31) 
(32) 

(33) 



(34) 
(35) 
(36) 

(37) 
(38) 
(39) 



(w). (40) 



The resulting determinant, despite its cumbersome expression, can be stud- 
ied by introducing the variable 



J ~ J„H' 



(41) 



which leads to 



D n (w) = 



48n(l-n 2 ) 
(n 2 - 4) 



!n( w )(y - yi)(y - 2/2X2/ - 2/3X2/ - 2/4), (42) 
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where 



n n n w n w 



(43) 




and hence 



(n 2 - 4) 



D n (w) = (i» + £/„h) (j» - £/„h) 



48n(l - n 2 ) 



■( 7 "H + (f + £)w)(/;H + (f-3^H). (44) 



3. First pair of scalar-mode ^-functions 

In our problem, the differential operator under investigation is the Lapla- 
cian on the Euclidean four-ball acting on metric perturbations. The bound- 
ary conditions for vector, tensor and ghost modes correspond to a familiar 
mixture of Dirichlet and Robin boundary conditions for which an integral 
representation of the ^-function and heat-kernel coefficients are immedi- 
ately obtained. New features arise instead from Eq. (44) , that gives rise to 
four different (- functions. On studying the first line of Eq. (44), we exploit 
the Cauchy integral formula to express the power — s of the eigenvalues and 
hence turn the £-function 



where 7 encloses the zeros X\,X2, ■■■,x oc of the function H n , which here 
equals J' n (x) ± ^J n (x). Such a combination of J„ and J' n is proportional to 
the power of degree ((3± — 1) of the independent variable multiplied by an 
infinite product, with (3 + {n) = n, /3_(n) = n + 2. Only the infinite product 
encodes information on the countable infinity of non-vanishing zeros, and 
hence one should divide xJ' n {x) ± nJ n (x) by x l3± . Last, rotation of contour 
to the imaginary axis (Dowker and Kirsten [22], Bordag et al. [27]), which 
brings in modified Bessel functions /„, jointly with setting w — zn, leads 
to the following integral formula: 



00 




into an integral, i.e. we use 





(sin7rs) 




'OO 



dzz 2s T-\og 



( 



znl' n {zn) ± nl n (zn) 



) 



7T 



z P±(n) 



(45) 
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The uniform asymptotic expansion of modified Bessel functions and 
their first derivatives (see Appendix) can be used to find (hereafter r = 
r(z) = (l + z 2 )-?) 

znUzn) ± nl n (zn) ~ ^=^(1 ± r) (l + g , (46) 

where (see Eqs. (139) and (141) in the Appendix for the functions Uk and 

Vk) 

p fc ,±(T) = (1 ± r)- 1 (v k {T) ± ™ fc (r)) , (47) 

for all fc > 1, and 

V fe=i / fe=i 

Thus, the ^-functions (45) obtain, from the first pair of round brackets in 
Eq. (46), the contributions (cf. Rcf. 22) 

A+( a ) = g„-(2.-2)(™) r^^^iogfi + ^ + ^-iy (49) 

A_ W S £ f d* log f 1 , 

^ 3 * •><> 9z ^ z 2 ; 

(50) 

where z 2 in the denominator of the argument of the log arises, in Eq. (50), 
from the extra z~ 2 in the prefactor in the definition (45). Moreover, 

the second pair of round brackets in Eq. (46) contributes X^iiA7,±( s )> 
having defined 

A, ± (s) = gn-(W)^£) r dz z -*°T jt± (T{z)) t (51) 
^3 * J ° dz 

where, from the formulae 

ri,±=pi,±, (52) 

T2,±=V2.,±-\vl±, (53) 

T 3 .± = P3,± - Pl,±P2,± + ^Pl,±i ( 54 ) 
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wc find 



^ = -iV 2 ^ 3+ r 4 ^ 5+ iV 6 < ^ 



3l± 128 8 + 640 T 12 + 128 8 1152 ' 1 j 
and hence, in general, 

3j 

r„ ± (r)=^.f^ ± )r a . (58) 

We therefore find, from the first line of Eq. (44), contributions to the 
generalized £- function, from terms in round brackets in Eq. (46), equal to 

oo 

X ±(s) =co (s)F±(s) + '£"j(s)F±(s), (59) 

where, for all A = 0,j (£r and being the Riemann and Hurwitz £- 
functions, respectively), 



u> x (s) = ^n-< 2s+A - 2 ) = Ch(2s + A-2;3) = 0?(2 S + A-2)-l-2-( 2s + A - 2 \ 

(60) 



n— 3 

while, from Eqs. (49)-(51), 



(61) 



#00 = ^ dz log(l + (1 + , 2 )-i) , 

F -( S) ^ _ 2 ^[?) jf"^ _ F + (S ) = -1-F+( S ), (62) 

^(-) S ^^i;L ± (-,a > 0)/W.±) > (63) 



<*=.J 

having set (this general definition will prove useful later, and arises from a 
more general case, where r a is divided by the &-th power of (1 ± r) in Eq. 
(58)) 

L^s.a.ft) = J T 2s+a (l-T)- s (l+T)- s (±b(l±T)- b - 1 -aT- 1 (l±T)- b yT. 

(64) 
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Moreover, on considering 



(65) 



and changing variable from z to r therein, all L-type integrals above can 
be obtained from 

Q(a,(3, 7 )= f 1 r a (l-Tf(l+rydr. 
Jo 

In particular, we will need 

L+(s) = -Q(2s,-s,-s-l), 



(66) 



(67) 



L + (s, a, b) = bQ(2s + a, -s, -s - b - 1) - aQ(2s + a - 1, -s, -s - b), (68) 

where, from the integral representation of the hypergeometric function, one 
has (Gradshteyn and Ryzhik [28]) 



Q{a > & 7) = r( r(o + y+2) 1) jF( ~ T ' « + + + 2; -i). 

For example, explicitly, 

^(.) = - r P'^-) r (. + l,2. + l,. + 2 ; -l). 
Now we exploit Eqs. (45), (46) and (59) to write 



(69) 



(70) 



[smns) 



J2n- (2s - 2) / dz 



n—3 



z-( 2s -V 
2(1 + z 2 ) 



+ nz- (2s+ ^ (VT + ^-l) 



Hence we find 



C+(0) = hm 



The first limit in Eq. (72) is immediately obtained by noting that 

limL+( S ) = -log(2), 

and hence 



(71) 



(72) 



(73) 



lim iuj {s)F+{s) = lim 

s— >0 s^O 



Ch(2 S -2;3)^^ 

7T 



0. (74) 
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To evaluate the second limit in Eq. (72), we use 

lim L + (s, a, 0) = — 1, 



(75) 



and bear in mind that u)j(s) is a meromorphic function with first-order 
pole, as s — > 0, only at j = 3 by virtue of the limit 



lim 



Cr(v) - 



Hence we find (see coefficients in Eq. (57)) 



OO / • \ oo 

X , sm7rs) x 

lim^ Wj ( S )f;+( S ) = Jim^i^^W 



= 7- 



(76) 



3j 



5> + ( S ,a,0)./^+) 



1 9 1 
0A^ Ja 79 



a— 3 



720' 



(77) 



while, from Eqs. (71) and (69), 



Hm(CUs) - X*00) = lim - 2; 3) + - 3; 3)) 



5 1079 

~ ~4 + 240 ' 

We therefore find, with the same algorithms as in Ref. 27, 

9 

' f(3, + ) _ 

45 ' 



(78) 



5 1079 1^/.:.. ^> 



4 240 2 



(79) 



a=3 



G(») 



240 



a— 3 



757 
~90~' 



(80) 



These results have been double-checked by using also the powerful analytic 
technique in Ref. 25. 



4. Further spectral asymptotics: elliptic and non-elliptic 
parts 

As a next step, the second line of Eq. (44) suggests considering C-functions 
having the integral representation (using again the Cauchy theorem and 
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rotation of contour as in Eq. (45)) 

oo 

E 



C | (s) = (™ff) v „ — 



7T 



n— 3 



9 

z -25 — log 

£72 



(znl' n {zn) + 



± n I n (zn) 



(81) 



To begin, we exploit again the uniform asymptotic expansion of modified 
Bessel functions and their first derivatives to find (cf. Eq. (46)) 



znl' n (zn)+ 



z 2 n 2 



± n I n (zn) 



n 2 e™ ? ' /l 



2^/2^m sfr \t 



r k,±(T) 



k=l 



where we have (bearing in mind that uo = vo = 1) 

2r / \ 

^,±(t) = Ufc(r) + ^— (( U fc-l( T ) ±TMfe-l(T)J, 

for all fc > 1. Hereafter we set 

ti n 

and rely upon the formula 



log(l + Q)~£ ( -l)^ 



k=i 



to evaluate the uniform asymptotic expansion (cf. Eq. (48)) 

fk,± (t(z)) \ ^ R ki ± (t(z)) 



log h+E 



fe=i 



E 

fe=i 



(82) 
(83) 

(84) 
(85) 
(86) 



The formulae yielding Rk,± from r k ,± are exactly as in Eqs. (52)-(54), with 
T replaced by R and p replaced by r (see, however, comments below Eq. 
(90)). Hence we find, bearing in mind Eq. (83), 



Ri ± = (1 T r)- 1 ( — t T -r 2 - -r 3 ± -r 4 
y 1 °, 8 24 24 



(87) 



" ; [ 16 8 16 4 16 8 16 ' ' 
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,_ 3 / 1721 3 441 4 597 5 1033 6 239 7 

■ v 7 V 384 128 320 960 80 

28 2431 221 _ 1105 1105 \ 
T 5 + 576 192 T 384 1152 J' 

and therefore 

fl ii± (r(z)) = (lTT)-^^V, (90) 

where, unlike what happens for the Tj t ± polynomials, the exponent of (lT T ) 
never vanishes. Note that, at r = 1 (i.e. z — 0), our rfe ! +(r) and Rk i+ (r) 
are singular. Such a behaviour is not seen for any of the strongly elliptic 
boundary- value problems [8] . This technical difficulty motivates our efforts 
below and is interpreted by us as a clear indication of the lack of strong 
ellipticity proved, on general ground, in Ref. 18. 

The Cb( s ) function is more easily dealt with. It indeed receives contri- 
butions from terms in round brackets in Eq. (82) equal to (cf. Eq. (50) 
and bear in mind that /3_ — j3 + = 2 in Eq. (81)) 

*_ W , J>- (2s - 2) jf * io g (&^) 



.(sin™) r,__ 2s d 1 _ 1, 



= w (s)- -/ dz z 2s — log^== = --wo(a), (91) 

7T Jo ^ v 1 + 2; 2 

and X^Jli Bj,-{ s )i having defined, with A = 0, j (cf. Eq. (51)) 

oo 

ux(s) = n-( 2s+x -V = ( H (2s + A - 2; 3), (92) 

n— 3 



, . (sin7rs) f 00 , 9c <9 . , 



(93) 



On using the same method as in Sec. 3, the formulae (81)-(93) lead to 



12 



„ , . 5 1079 5 1 A > 206 

c ^-4 + w + ri6^ c «'^¥' (94) 

a— 3 

a result which agrees with a derivation of (0) relying upon the method 
of Ref. 25. 

Although we have stressed after Eq. (90) the problems with the Cb( s ) 
part, for the moment let us proceed formally in the same way as above. 
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Thus wc define, in analogy to Eq. (91), 



.(sin™) r,_- 2s d^j 1 



B + ( S )^. 0{S )^-^J dzz--lo S ^-r(z)), (95) 

and, in analogy to Eq. (93), 

B^^.^L—Jjf dzz-^-R^ + (r(z)). (96) 

In order to make the presentation as transparent as possible, we write out 
the derivatives of Rj,+- On changing integration variable from z to r we 
define 

Cj{r) = ^RjAt), (97) 
and we find the following results: 

C 2 (r) = (1 - ,)"" ( - %r + ^ - ? T 3 + ^ _ ^ _ 33 6 + « r 



^r 8 ), (99) 



/1721 2 
( 128 T 


441 , 


1635 4 

T — 


163 5 




~ *T + 


128 


TT 


221 9 


5083 in 


1105 n 


1105 


16 T 


T + 

128 


32 T 


128 



4223 8 _ ^ g _ uuo. 10 n _ ^ 12 | ( } 

64 16 128 32 128 j' v ; 

so that the general expression of Cj (r) reads as 

4j 

Q(r) = (1 - r)-^ 1 J] X«r", Vj = 1, . . .,oo . (101) 

a=j-l 

These formulae engender a Cb(0) which can be defined, after change of 
variable from z to t, by splitting the integral with respect to r, in the 
integral representation of Cs( s )> according to the identity 

1 rjj, /-l 

dr — dr + I dr, 
o Jo J n 

and taking the limit as /z — > 1 a/ter having evaluated the integral. More 
precisely, since the integral on the left-hand side is independent of fi, we 
can choose fi small on the right-hand side so that, in the interval [0, /i] (and 
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only there!), we can use the uniform asymptotic expansion of the integrand 
where the negative powers of (1 — r) are harmless. Moreover, independence 
of /i also implies that, after having evaluated the integrals on the right- 
hand side, we can take the \i — > 1 limit. Within this framework, the limit 
as fi — ► 1 of the second integral on the right-hand side yields vanishing 
contribution to the asymptotic expansion of Cb( s )- 
With this caveat, on defining (cf. (66)) 



Q„(a,/J, 7 ) = / T a (l-Tf(l + rydT, (102) 
Jo 

we obtain the representations 
f sin it s ) " 

B+(s) = -wo(*)- L - Q A( (2s, -s - 1, -s) + Q^ls, -s, -s - 1) 

7T L 

-Q„(2«-l, -a,-*)], (103) 



Bj, + (s) = J] ^^(2^ + «, -« - i - 1, s). (104) 

o=J — 1 

The relevant properties of Q jU (a,/3, 7) can be obtained by observing that 
this function is nothing but a hypergeometric function of two variables [28] , 
i.e. 

u a+1 

Q»{a,f3,j) = ^—F 1 (a + l,-(3,- 1 ,a + 2- l i,- l i). (105) 
a + 1 

In detail, a summary of results needed to consider the limiting behaviour 
of Cb( s ) as s is 

, , (sin7rs) „, , 

w (s)- -~-5s + 0(s 2 ), (106) 



w .( s) l!^ „ 1^. 3 + b jA s + 0(s 2 ), (107) 

7T Z 

lim Q M (2s, -s - 1, -s) ~ -- + O(s ), (108) 



lim Q (2s, -s, -a - 1) ~ log(2) + O(s), (109) 



lim Q^(2s - 1, -s, -a) ~ + O(s), (110) 

(ti— »1 zs 
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lim QJ2s + a, -s - j - 1, -s) 
r(-j-s)r(a + 2s+l) 



2 Fi (a + 2s + 1, s, a - j + s + 1; -1) 
+ 6j,o(a) + 0(«), (111) 



T(a -j + s + 1) 
bj.,-i(a) 



where 



hi = "I - 22 ^' + CflO' - 2)(1 - ^,3) + 7^,3, 
(-1)^ +1 r(a+l) 



&j,_i(a) = 



j! r(a-j + l) 



jt(1 - 5a,j-l), 



(112) 
(113) 



and we only strictly need 63.0(a) which, unlike the elliptic cases studied 
earlier, now depends explicitly on a and is given by (ip being the standard 
notation for the logarithmic derivative of the T-function) 



. , . ir( + i) 

t»3,o(a) = c 



6T(o-2) 
- V(a-2)-V(4) 



- log(2) - i(6a 2 - 9a + l)^-fy + 2</>(a + 1) 



(114) 



Remarkably, the coefficient of ^ in the small-s behaviour of the gener- 
alized function Cb( s ) is zer0 because it is equal to 
12 1 12 

lim S <+(*) = ^6 3 ,-i(a)^ = -5>(a-l)(a-2)X( 3 ), (115) 

a=2 a=3 

which vanishes by virtue of the rather peculiar general property 

E r^t+n ^ = E fiV - 0^ - 0, v, = 1, . . ., 00, (lie) 

a=j ^ ■' ' a=j 1=0 

and hence we find eventually 



£(0) - 



5 

4 + 



1079 


5 


1 12 

a=2 


+ 240 


+ 2 


1079 
240 + 


599 


296 


720 


" 45"' 



00 4j 

E^ E ' 

i=l a=j'-l 



(117) 



because the infinite sum on the first line of Eq. (117) vanishes by virtue 
of Eqs. (113) and (116), and exact cancellation of log(2) terms is found to 
occur by virtue of Eq. (116). 
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To cross-check our analysis, we use Eq. (83) to evaluate 

4r 2 



and hence we find 



64 

+ 



^2,+ — R2, 
4r 2 



rk,+ (r) - r- fe ,-(r) 
Ri,+ = Ri, 

4r 2 



(1-r 2 ) 
4r 2 



Wfe-i(r), 



(1 



(1-r 2 ) 



(1-r 2 ) 
4r 2 



2r 2 



Ui 



"2 



2r 2 



(1-r 2 ) 



(1-r 2 ) 

Ml — M1-R1 _ 



Ri, 



R2 



4r 2 



(1-r 2 ) 



(118) 

(119) 
(120) 
Ri.- 



3 (1-r 2 ) 3 (1-r 2 ) 

and so on. This makes it possible to evaluate Bj_ + (s) — Bj_(s) for all 
j = 1, 2, ...00. Only j = 3 contributes to Cs (0) (see below) and we find 



(121) 



B3,+(*)--B 3 ,_(«) = -ws(«) 



'sin 7rs) 



• lim rd T r 2s (l-r)- s (l+r)- s |-(i?3, + -ii3,-) 



(122) 



80r 3 - 24r 5 + 32r 7 - 8r 9 , (123) 



The derivative in the integrand on the right-hand side of Eq. (122) reads 
as 

^( J R 3 ,+ - J R3,-) = (l-r)- 4 (l + r)- 4 

and hence we can use again the definition (102) and the formula (105) 
to express (122) through the functions Q M (2s + a, — s — 4, — s — 4), with 
a = 3,5, 7, 9. This leads to 

£(o) = Ci(o) + i?3, + (o)-i?3,-(o) 
4 



- rf0) Ly- £([+!) 

24f^T(i-2) 



V(Z + 2) 



1 



.(3) 
V 2Z+1 



206 „ 296 
h 2 = , 

45 45 



(124) 



where >4z+i are tne f° ur coefficients on the right-hand side of (123). Regu- 
larity of Cb( s ) & t the origin is guaranteed because lim s ^ s Cb( s ) ^ s propor- 
tional to 

1 r(i + i) {3 ) 



K 2Z+1 — 0' 
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which is a particular case of the peculiar spectral cancellation (cf. (116)) 



ama*(j) r ( {a+1) ) 



a— a n 



a(j) 1 V 2 



where a takes both odd and even values. The case j = 3 is simpler because 
then only Ka ^ coefficients with odd a are non-vanishing. 

Remaining contributions to C(0), being obtained from strongly elliptic 
sectors of the boundary-value problem, are easily found to agree with the 
results in Rcf. 19, i.e. 

278 

£(0) [transverse traceless modes] = — (126) 

494 

C(0) [coupled vector modes] = (127) 

15 

C(0) [decoupled vector mode] = — — , (128) 

C(0) [scalar modes(ai, e\\ a,2, b 2 , e 2 )] = —17, (129) 

149 

((0) [scalar ghost modes] = — (130) 

77 

C(0)[vector ghost modes] = — , (131) 

5 

C(0) [decoupled ghost mode] = -. (132) 
Our full ((0) is therefore, from (79), (80), (94), (117), (126)-(132), C(0) = 

142 
45 ■ 



5. Concluding remarks 

We have studied the analytically continued eigenvalue conditions for met- 
ric perturbations on the Euclidean four-ball, in the presence of boundary 
conditions completely invariant under infinitesimal diffcomorphisms in the 
de Donder gauge and with the a parameter set to 1 in Eq. (10). This 
has made it possible to prove that only one sector of the scalar-mode de- 
terminant is responsible for lack of strong cllipticity of the boundary-value 
problem (see second line of Eq. (44) and the analysis in Sees. 3 and 4). 
The first novelty with respect to the work in Ref. 18 is a clear separation 
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of the elliptic and non-elliptic sectors of spectral asymptotics for Euclidean 
quantum gravity. We have also shown that one can indeed obtain a regular 
C-function asymptotics at small s in the non-elliptic case by virtue of the 
remarkable identity (116). Our prescription for the £(0) value differs from 
the result first obtained in Ref. 19, where, however, neither the strong el- 
lipticity issue [18] nor the non-standard spectral asymptotics of our Sec. 4 
had been considered. 

As far as we can see, the issues raised by our results are as follows. 

(i) The integral representation (81) is legitimate because the second line of 
Eq. (44) corresponds to the eigenvalue conditions, for n > 3, 

F£{n,x) = J' n {x) + (-| ± J n (x) = 0. (133) 

For both choices of sign in front of ^, if Xi is a root, then so is — Xi, 
with positive eigenvalue Ei = xf (having set the 3-sphere radius q = 1 for 
simplicity). For any fixed n, there is a countable infinity of roots Xi and they 
grow approximately linearly with the integer i counting such roots. The 
function Fg admits therefore a canonical-product representation (Ahlfors 
[29]) which ensures that the integral representation (81) reproduces the 
standard definition of generalized ^-function, i.e. 

C(s) = £ d{E k )E-% 

E k >0 

where d(Ek) is the degeneracy of the eigenvalue E^. 

(ii) Even though the lack of strong ellipticity implies that the functional 
trace of the heat semigroup no longer exists, and hence the Mellin transform 
relating ^-function to integrated heat kernel cannot be exploited, it remains 
possible to define the functional determinant of the operator P acting on 
metric perturbations. For this purpose, a weaker assumption provides a 
sufficient condition, i.e. the existence of a sector in the complex plane free 
of eigenvalues of the leading symbol of P (Secley [30]). Note also that, if 
one looks at the A\ heat-kernel coefficient for boundary conditions involv- 
ing tangential derivatives [8], it is exactly for the ball that the potentially 
divergent pieces involving the extrinsic curvature in A\ cancel. Thus, on 
the Euclidean ball cancellations take place that maybe could explain why 
C(0) is finite. This might be therefore a very particular result for the ball. 

(iii) By virtue of standard recurrence relations among Bessel functions, 
the eigenvalue conditions (133) are equivalent to studying the eigenvalue 
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conditions 



Fg(n,x) = J n {x) =F -J„-i{x) = 0, 



(134) 



where the eigenvalues E(i,n,±) are obtained by squaring up the roots 
x(i,n, ±). The equation for Fg(n 7 x) can be further re-expressed in the 
form 



The functions F B differ therefore by one term only and this term gets 
small as x gets larger. The numerical analysis confirms indeed that a p(i, n) 
positive and much smaller than 1 exists such that one can write (Esposito 
et al. [31]) 

E(i, n, +) = E(l, n, +)5 hl + E(i - 1, n, -)(1 + p(i, n))(l - 6 iA ), (136) 
for all n > 3 and for all i > 1. 

(iv) The remarkable factorization of eigenvalue conditions, with resulting 
isolation of elliptic part of spectral asymptotics (transverse-traceless, vector 
and ghost modes, all modes in finite-dimensional sub-spaces and three of 
the four equations for scalar modes), suggests trying to re- assess functional 
integrals on manifolds with boundary, with the hope of being able to obtain 
unique results from the non-elliptic contribution. If this cannot be achieved, 
the two alternatives below should be considered again. 

(v) Luckock boundary conditions (Luckock [32]), which engender BRST- 
invariant amplitudes but are not diffcomorphism invariant [15]. They have 
already been applied by Moss and Poletti [12, 33]. 

(vi) Non-local boundary conditions that lead to surface states in quantum 
cosmology and pseudo-differential operators on metric and ghost modes 
(Marachevsky and Vassilevich, Esposito [34]). Surface states are particu- 
larly interesting since they describe a transition from quantum to classical 
regime in cosmology entirely ruled by the strong ellipticity requirement, 
while pseudo-differential operators are a source of technical complications. 

There is therefore encouraging evidence in favour of Euclidean quantum 
gravity being able to drive further developments in quantum field theory, 
quantum cosmology and spectral asymptotics (see early mathematical pa- 
pers by Grubb [35], Gilkey and Smith [36]) in the years to come. 




(135) 
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Appendix: Olver expansions 



In Sees. 3 and 4 we use the uniform asymptotic expansion of modified 
Bessel functions I v first found by Olver [37]: 



I v {zv) 



r i + E 



Uk(r) 



k=l 



where 



T = (l + Z 2 )-2, 77= (1+Z 2 )* +l0g 



vT 



(137) 



(138) 



This holds for v — > oo at fixed z. The polynomials Ufe(r) can be found from 
the recurrence relation [27] 



«Hi(T) = ^ 2 (l-r 2 K(r) + i 



dp (1 - 5p 2 )u k (p), 



(139) 



starting with Uo(t) — 1. Moreover, the first derivative of has the follow- 
ing uniform asymptotic expansion at large v and fixed z: 



e"" (l + z 2 )i 

J2-KV z 



Vfe(r) 



fc=i 



with the i;^ polynomials determined from the Mfc according to [27] 



Ufe(r) = Ufe(r) + r(r 2 - 1) 
starting with vq{t) = Uq(t) = 1. 



-Ufc-i(r) +TM / fe _ 1 (r) 



(140) 



(141) 
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